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Abstract 

The influence of nuclear polarization on the bound-electron g factor in heavy 

hydrogenlike ions is investigated. Numerical calculations are performed for 

the K- and L-shell electrons taking into account the dominant virtual nuclear 

excitations. This determines the ultimate limit for tests of QED utilizing 

measurements of the bound-electron g factor in highly charged ions. 
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Recent high-precision experiments for measuring the bound-electron g factor in hydro- 
genlike carbon have reached a level of accuracy of about 2 x 10 -9 []I]|§. As a consequence, 
this has led to a new independent determination of the electron mass H . Via investigations 
of the g factor of a bound electron in a highly charged ion one can probe nontrivial effects 
in bound-state QED as sensitive as in high-precision Lamb shift experiments. A further 
improvement in accuracy and the extension to systems with higher nuclear charge numbers 
Z up to hydrogenlike uranium is intended in the near future jy. Studies of g factors in 
heavy ions are of particular importance, since they can provide a possibility for an indepen- 
dent determination of the fine-structure constant @||, nuclear magnetic moments ||, and 
nuclear charge radii. In order to achieve a level of utmost precision in corresponding theo- 
retical calculations, one has to account for the relativistic, higher-order QED, nuclear-size, 
nuclear-recoil, and nuclear-polarization corrections PHJ^I- Investigations of QED effects in 



heavy systems are strongly restricted by the uncertainty due to the finite nuclear size [§,14 



In Ref. a specific difference has been introduced for bound-electron g factors in H- and 
Li-like ions, for which the uncertainty due to the nuclear-size effect can be significantly re- 
duced. With an apparent accuracy of 10~ 9 for the bound-electron g factor one could probe 
higher-order QED corrections even for uranium ions, provided that nuclear polarization 
effects remain negligible. 

In the present Letter, we evaluate nuclear-polarization corrections to the g factor in 
hydrogenlike ions. This reduces the remaining source of uncertainties in the prediction of 
nuclear effects. Moreover, we determine the ultimate limit of accuracy for QED tests in 
measurements of the bound-electron g factor in highly charged ions. 

We consider a hydrogenlike ion with a spinless nucleus, which is placed in a homogeneous 
external magnetic field 7i corresponding to a vector potential A(r) — [H x r]/2. The energy 
shift of a bound-electron level a within first-order perturbation theory in the magnetic field 
(see Fig. |l|(a)) is given by (h — c — 1) 

AE a = (ij a \V H m, (1) 
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where 

K« = M(7f.[rxa]). (2) 



Choosing the z axis along the direction of the field 7i, i.e., 7i = (0,0, H), one obtains for 
the energy shift 

AE a = fioHmjga, (3) 

where /io = |e|/(2m) is the Bohr magneton, rrij is the ^-projection of the total angular mo- 
mentum, and g a is the bound-electron g factor, which depends on the electron configuration. 
In the case of a Dirac electron in the Coulomb field of an infinitely heavy point-like nucleus, 
it yields |L6 



A = JoTTj («£ -5) ■ w 

Here k = (j + 1/2) (— l)J+' +1 / 2 is the relativistic angular-momentum quantum number, j is 
the total angular momentum of the electron, I = j ± 1/2 defines the parity of the state, e nK 
is the one-electron energy of the state given by 

l + n r 

£uk = m — , (5) 

where n r = n — \k\ is the radial quantum number, n is the principal quantum number, 
7 = J k 2 — (aZ) 2 , N = J{n r + ^j) 2 + (aZ) 2 , and a = e 2 is the fine-structure constant. 



Due to various QED and nuclear effects, the observed bound-electron g factor deviates from 
its Dirac value (^). Here we consider the nuclear-polarization correction Ag nK , which is of 
particular importance in heavy ions. 

The dominant contribution to the nuclear-polarization effect for heavy nuclei arises from 
virtual collective nuclear excitations. Three types of collective modes should be taken into 
account: (a) rotations of the deformed nuclei; (b) harmonic surface vibrations; and (c) gi- 
ant resonances. In Ref. flI7| , a relativistic field theoretical approach has been developed, 
where the nuclear-polarization effects are treated perturbatively, incorporating the many- 
body theory for virtual nuclear excitations within bound-state QED for atomic electrons. 



The contribution of the nuclear vector current can be omitted, because the velocities asso- 
ciated with collective nuclear dynamics are nonrelativistic jTj|. Accordingly, one is left with 
the longitudinal component of the effective photon propagator Dqq only due to the nuclear 
transition density-correlation function. In Coulomb gauge, it can be represented in terms of 
a multipole decomposition as follows |T7 



D 00 (r, r'; co) = £ B{EL) ^ ^f^J Y^) ■ Y L (p,')). (6) 

Here ujl = El — Eq are the nuclear excitation energies with respect to the ground-state 
energy E Q of the nucleus and B(EL) = B(EL;0 — > L) are the corresponding reduced 
electric transition probabilities. The radial shape parametrizing the nuclear transitions is 
carried by the functions 

4tt 



FrXr) 



(2L + 1)R% [R L 



r L R% 



e(i2b-r) + ^e(r-i2o) 



(7) 



for the case of multipole excitations with L > 1 and 



F (r) 



Q(Ro-r) (8) 
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for monopole excitations, respectively. Here Ro is an average radius assigned to the nucleus 
in its ground state. The presence of O-functions in the expressions (^) and (§) reflects 
the sharp surface approximation for collective excitations. The form (§) of the propagator 
is convenient for numerical evaluations, since the parameters characterizing the nuclear 
dynamics ujl and B(EL) can be taken from experiment. Nuclear-polarization corrections 
to the Lamb shift (see graph on Fig. 0(b)) have been calculated in Refs. [p| -pO|. 

The nuclear-polarization contribution to the bound-electron g factor appears as the 
lowest-order nuclear-polarization correction to the diagram |l](a). To first order in V^, this 
perturbation gives rise to a modification of the wave function, of the binding energy, and of 
the electron propagator. The corresponding contributions are referred to as the irreducible 
part, the reducible part, and the vertex part, respectively. The nuclear-polarization en- 
ergy shift of the state under consideration may be represented by the diagrams depicted in 



Figs. |I](c) and 0(d). Let us consider first the energy correction due to the irreducible part 
of the graph |I|(c) 

AE? = -2a £ B(EL)^-'t T~ ^ F ^"> 
tm '2L + 1 ft J2me a -u-B n (l-i0) 



X 



{ii> n \F L Y£ M \ip k ) & k \V H M 



u 2 -u 2 L + iO e a -£ k (l-iO)' 
Here the indices n and k in the sum run over the entire Dirac spectrum. 

After integration over frequencies u and summation over angular projections, Eq. 
takes the form 



(9) 



(nn\FL\niK,i) (niK>i\FL\n2n) (ri2K\ra x \nK) 



(10) 



£-riK &n\n\ Sgn^£ niKl Jto>£ ^nn ^ri2H 

The sum over Ki is restricted to those intermediate states, where I + l\ + L is even. In 
Eq. (|T0"D, a two- component radial vector (r\nn) is determined by 

'Pn K (r)\ 



(r\nn) 



^Qnft(^) 



(11) 



where P nK (r) = rg nK (r) and Q nK (r) = rf nK (r), with sw(r) and / nre (r) being the upper and 
lower radial components of the Dirac wave function | 2l| , respectively. The radial matrix 
element is given by 



(a\F L \b) = J drF L (r) [P a (r)P b (r) + Q a (r)Q b (r)} 



(12) 



and a T denotes the Pauli matrix. The sum 



^ (r\n K)(n K\ra x \nK) 

{r\nK) = 2^ 



(13) 



can be evaluated analytically using the generalized virial relations for the Dirac equation 



22[. The upper and lower components P nK {r) and Q nK (r) of the vector {r\nn) read [JTO 
P n Jr) = ( inn — + KS nK ) r + olZk Q nK (r) + 
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(14) 
(15) 



Finally, the corresponding irreducible part of the correction Ag nK can be expressed as 

2 (nK,\FL\niKi) (niKi|Fi|?m) 



a Km 



(~iJl 2 



- sgn(e ni/s Jw L 



(16) 



The reducible part of the graph depicted in Fig. |l|(c) has to be considered together with 
the contributions resulting from diagrams [l|(a) and|l|(b). The corresponding corrections to 
the energy shift reads 

~° duj (Tp a \F L Y LM \ip n ) 
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leading to the ^-factor correction 
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(18) 



sgn^„ 1K ju; L j 

Here is the Dirac g factor given by Eq. (|J). In Eqs. QlfiD and fllSD , the sum Z + Zi + L 
again should be even. 

Let us now turn to the nuclear-polarization correction to the vertex as depicted in 
Fig. |l](d). The corresponding energy shift is determined by 



AEr 



ftr 21 + 1^ J 

L,M n,fc_„ 
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ijj- 



(19) 



uj\ + iOe a - uj - e k (l - iO)" 

The integration over and the summation over angular variables leads to the corresponding 
expression for Ag™ r , which is conveniently represented as the sum of a pole term 
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and of a residual term 



sgn(£ niK1 )o; L ] 5 



(20) 
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Vj(j + 1 )( 2 J + 1 ) L ni,n s «i,K 2 2 2 2 



x f Ji J2 1 1 (n/t|F L |ni/ti)(n 2 K 2 |F L |nK) (ni/ti|ro- z |n 2 ft:2) 

respectively. Here AgJ^ 1 accounts for the terms with rii = n 2 and Ki = k 2 in the sums over 
intermediate states. The prime in the sum in Eq. indicates that e niKl ^ s n2K2 when 
Ki = k 2 , i.e., the pole contribution is supposed to be omitted. In Eqs. ( p0|) and (pTf), the 
value I + h + L has to be even. A second condition in Eq. ( pi]) is that the sum li + l 2 should 
be even as well. The total nuclear-polarization contribution to the g factor is determined 



by the sum of all contributions given by Eqs. (|T6|), (|18|), (p0|), and (pT|). 

We have evaluated the nuclear-polarization correction to the g factor taking into account 
a finite set of dominant collective nuclear excitations (see Table |) . For low- lying rotational 
and vibrational levels, the corresponding nuclear parameters, u>l and B(EL), have been 
taken from experiments on nuclear Coulomb excitation. In our estimates for the contri- 
butions due to giant resonances, we utilized phenomenological energy-weighted sum rules 
UH. Th e latter are assumed to be concentrated in single resonant states. In the present 



calculations, contributions due to monopole, dipole, quadrupole, and octupole giant reso- 
nances have been taken into account. To evaluate the infinite summations over the entire 
Dirac spectrum, the finite basis set method has been employed. Basis functions have been 
generated via B splines including nuclear-size effects ||24)| . The major contribution to the 
Ag nK results from the correction to the wave function (|T6|). This is due to the fact that the 
matrix element of the atomic magnetic-moment operator is saturated over the scale of atomic 
distances, while the influence of nuclear-polarization is essential in the vicinity of nucleus 
only. In the irreducible term, atomic and nuclear scales come into play simultaneously. 

According to our numerical results, we conclude that introducing a specific difference of 
g factors of H- and Li-like heavy ions |l5j does not eliminate all the nuclear effects. The 
influence of intrinsic nuclear dynamics becomes noticable at a level of accuracy of about 
10~ 9 for nuclei in the medium Z-range and increases up to 10~ 6 in uranium. Since nuclear- 
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polarization effects set a natural limit up to which bound-state QED can be tested, one is 
faced here with a situation similar to the one in Lamb shift experiments. However, within 
the expected accuracy of 10~ 9 in g-factor experiments with heavy highly charged ions one 
may provide a tool for probing internal nuclear structure and for testing specific nuclear 
models. 
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grateful for financial support from RFBR (Grant No. 01-02-17246) and from the Alexander 
von Humboldt Foundation. G.P. and G.S. acknowledge financial support from BMBF, DFG, 
and GSI. 
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FIGURES 



(a) (b) (c) (d) 

FIG. 1. Diagrams representing the interaction of a bound electron with the external magnetic 
field (a), the lowest-order nuclear-polarization effect (b), and the nuclear-polarization correction 
to the bound-electron g factor (c) and (d). The heavy line denotes the nucleus. The contribution 
corresponding to the graph (c) should be counted twice. 



l~V/\/\/V<| 
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TABLES 

TABLE I. Nuclear-polarization effects to the g factor of K- and L-shell electrons in hydrogenlikc 
ions. Column (a): contributions from low- lying rotational and vibrational nuclear modes using 
experimental values for nuclear excitation energies u>l and electric transition strengths B(EL); (b) 
contributions from giant resonances employing empirical sum rules [20,23]; (c) total effect. The 
numbers in parentheses are powers of ten. 
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